QUANTUM A r Q-SYSTEM SOLUTIONS AS Q-MULTINOMIAL SERIES 



PHILIPPE DI FRANCESCO 

Abstract. We derive explicit expressions for the generating series of the fundamental so- 
lutions of the A r quantum Q-system of Ref. [P. Di Francesco and R. Kedem, Noncommuta- 
tive integrability, paths and quasi- determinants, preprint arXiv : 1006 . 4774 [math-ph] ] , 
expressed in terms of any admissible initial data. These involve products of quantum 
multinomial coefficients, coded by the initial data structure. 
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I. Introduction 

The study of some discrete integrable systems, taking the form of recursion relations, 
i.e. evolution equations in a discrete time n e Z with suitable conservation laws, has 
recently shed some new light [3] on the positivity conjecture of cluster algebras [8j. Indeed, 
admissible sets of initial data for such systems are particular clusters in some specific, non- 
finite type cluster algebras, whereas cluster mutations are implemented by the update of 
initial data via local application of the evolution equation. The positivity conjecture for 
cluster algebras then boils down to the following property: the solutions of such systems 
are Laurent polynomials with non-negative integer coefficients of any of their admissible 
initial data. 
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In [3], this was proved for the so-called Q-system for A r , by expressing solutions as 
partition functions for weighted paths on some target graphs, the weights being explicit 
Laurent monomials of initial data. An extremely useful tool for such representations is the 
notion of (multiply branching) continued fraction. We were able to prove that mutations of 
initial data are implemented by local rearrangement of the continued fraction expressions 
for the generating series of fundamental solutions of the Q-system. This was then extended 
to the so-called T-systems in (5] for initial data forming periodic stepped surfaces, and 
finally to the most general initial data in [2] , by use of a manifestly positive network path 
formulation. In all these cases, positivity follows from a form of discrete path integral 
representation of the solutions. 

The cluster algebra structure has a natural quantum version pQ , in which cluster variables 
obey quantum commutation relations within each cluster. This led to the natural definition 
of quantum Q-systems [TJ. In an analogous spirit, it was shown in j5] that the T- system may 
be viewed as a Q-system, involving non-commutative (time-ordered) variables. An even 
broader non-commutative version is known for the fully non-commutative A\ Q-system, for 
which Laurent positivity was conjectured by M. Kontsevich, and subsequently proved in 
[5] . The latter proof relies on an extension of the previous path formulations to paths with 
non-commutative step weights: the partition function of such paths is the sum over paths 
of the product of step weights taken in the same order as the steps are taken. In [TJ, such 
paths were used to investigate non-commutative versions of the Q-systems. In particular, 
a compact formulation in terms of non-commutative continued fractions was obtained. 

Except in very specific cases, very few explicit expressions of cluster variables in terms 
of fundamental data are known. For the (classical) A r Q-system, such expressions were 
derived in jl] for the generating series of its fundamental solutions. 

The aim of this note is to generalize these expressions for the generating series of the 
fundamental solutions of the quantum A r Q-system, by using the non- commutative con- 
tinued fraction expressions of [7J. The results are summarized in our main Theorem 13.121 
below, which expresses these generating functions for any admissible initial data as explicit 
series with coefficients that are Laurent polynomials with coefficients in Z + [g, where 
q is the parameter of the quantum deformation. 

Acknowledgments: We thank R. Kedem for helpful discussions and the Mathematical 
Sciences Research Institute, Berkeley, for hospitality during the program "Random Matrix 
Theory, Interacting Particle Systems and Integrable Systems" (fall 2010) during which this 
work was initiated. 

2. The quantum A r Q-system: definitions 

2.1. The system. Let A« = i(r + 1 — i), i = l,2,...,r and q G C*. The A r quantum 
Q-system [TJ is an evolution equation for variables Rij, « G [l,r] and j G Z, elements of a 
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non- commuting unital algebra: 

(2.1) g A * ^j+i^j-i = Rlj + Ri+ijRi-ij (i G [1, r],j G Z) 

with Rqj = R r +ij = 1 for all j G Z. 

2.2. Initial data. This is a three-term recursion relation in the variable j, which allows to 
determine all Ri j in terms of any initial data covering two consecutive values of j. Initial 
data are indexed by Motzkin paths m = (m 1 , ...,m r ) with m i+1 -m, G {0, 1, —1}. They read 
x m = (Ri,mi, Ri,mi+i) cn n an d are transformed into each-other via (forward/backward) 
mutations /i^ 1 that act on the Motzkin paths via (//f(m)) = rrij + eSj^, e — ±, whenever 
the result is itself a Motzkin path. The fundamental initial data corresponds to the null 
Motzkin path m = (0, 0, 0). 

Within each such set of initial data, the variables R obey the following commutation 
relations: 

(2-2) Ri,jRk,m ^ *'* Rk,mRi,j 

where 

(2.3) A i>k = Min(i, k) (r + 1 - Max(i, fc)) 

2.3. Commuting limit. Setting q — 1, in which case all variables commute, we 
recover the (commuting) A r Q-system: 

(2.4) R^R^ = R\ 3 + />',. ,.,/.', (^ G [1, r], j G Z) 
with _R ,j — R r +ij = 1 f° r ah j £ Z. 

2.4. Quantum cluster algebra for the A r Q-system. 

2.4.1. Cluster algebra and quantum cluster algebra. A cluster algebra of finite rank n with- 
out coefficients is a commuting algebra generated by invertible variables forming n-vectors 
~x.it) = (xi(t),X2(t), ...,x n (t)) attached to the vertices t of an infinite n-valent tree with edges 
labeled 1,2, ...,n around each vertex. Vectors x(t) and x(tt) corresponding to vertices t,u 
connected by an edge labeled k are related by a mutation relation of the form: 

Xi(u) = Xi{u) for i 7^ k 

n n 

(2.5) x k (u) Xk (t) = n^w [Bi ' feWi+ +n^w hBi ' feW]+ 

i=\ i=l 

where B{t) is an n x n skew-symmetrizable matrix, called the exchange matrix, with entries 
in Z, attached to the vertex t, and subject to the mutation relation 

Bij(u) = —Bij(t) if % = korj = k 

(2.6) Bij(u) = B hj {t) + sgn{B hk {t))[B hk {t)B k ^t)] + otherwise 
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where [x] + = Max(x, 0). In the following we'll be dealing only with cluster algebras with 
skew-symmetric exchange matrices. In that case, we may represent each matrix B(t) as a 
quiver with n vertices corresponding to the cluster variables, and with B(t)ij = the number 
of arrows from vertex i to vertex j whenever B(t)ij > 0. 

A cluster algebra is entirely specified by the pair (x(t ), B(t )) of cluster variables and 
exchange matrix at an initial vertex to, also called fundamental seed. 

The cluster algebras have the Laurent property, that any cluster variable Xk{t) may be 
expressed as a Laurent polynomial of the cluster variables at any other vertex u of the tree. 
It was conjectured in [5] that these polynomials have non-negative integer coefficients. 

A quantum cluster algebra [1] of rank n is a non-commuting version of the former de- 
fined as follows. Starting from an ordinary cluster algebra data, we introduce an extra 
n x n integer matrix A(to), forming with -B(to) a "compatible pair", namely such that 
A(to)B(to) = d, a diagonal matrix with positive integer entries. The matrix A(to) en- 
codes the quantum commutation relations obeyed by the initial cluster variables x(to), 
with Xi(t )xj(t ) = g Ai '^'°^Xj(t )^i(^o), where q is a fixed central element of the algebra. 
The mutation of cluster variables is defined via an analogous, non-commuting formula, 
while that of exchange matrices remains the same (12.61) . Compatibility fixes A(t) for every 
vertex t as well. 

Quantum cluster algebras also satisfy an analogous Laurent property. The positivity 
conjecture claims that the coefficients of the Laurent polynomials belong to g -1 ]. 



2.4.2. Cluster algebra for the commuting A r Q-system. The cluster algebra for the (com- 
muting) A r Q-system [9] has rank 2r, and a fundamental seed made of the cluster xo = 
x m — (Rifli -^2,0, Rr,o, Ri,i, Rr,i), an d of the 2r x 2r skew-symmetric exchange matrix 

Bq = B ma — ( ), C the Cartan matrix of A r , with entries 



(2.7) = 2^- - (i,je[l,r}) 

All the initial data x m are clusters in this cluster algebra. They are obtained from xo via 
iterated cluster (forward or backward) mutations of the form [if above, namely leaving all 
cluster variables unchanged except R i>mi Ri,m l+ 2 = (Ri, mi +i +%,m i +i^-i,m i +i)/^,mi 
(/i+) or Ri >m+1 -»■ Ri, mi -i = {Rl mi + Ri+x, mi Ri-i,mi) I 'Ri, mi +i (A r )> when all three terms 
are cluster variables in the original cluster. With our choice of fundamental seed, the first 
r variables always have even indices m while the r next have odd ones. The following is a 
sequence of forward mutations applied successively on the fundamental initial data x mo in 
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the case r = 3: 
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Here, each dot (m, a) G Z 2 represents a cluster variable R a , m , and the corresponding 
Motzkin path is the set of leftmost dots in the pairs. 

The (skew-symmetric) exchange matrix B m for each Motzkin path m was computed 
explicitly in [3]. As explained above, it can be represented as a quiver with 2r vertices 
indexed by the initial data indices (rrii,i) and (m« + l,z). Here's the example for r = 3: 



Ar 



7° 

to 





KIN 



where we have represented the vertices on the same Z 2 grid (m, a) as for initial data above, 
and where arrows indicate various forward mutations of the corresponding index. Here, 
we have only represented the exchange matrices for a fundamental set of Motzkin paths 
modulo a global translation by (1,0), namely the set M r = {(m a ) ae [i ir ]|Min a (m a ) = 0}. 
The exchange matrix is actually quasi-periodic, namely: B m+ i = —B m , where we denote 
by m + 1 the Motzkin path (mi + 1, 1712 + 1, m r + 1). 

More generally, the exchange matrix B m is constructed as follow^ Given the Motzkin 
path m = (m a ), we simply represent its vertices and their translates by the vector (1, 0) on 
the (m, a) plane, and we represent either of the three following local arrow configurations, 



This construction is due to R. Kedem. 
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depending on whether the Motzkin path is locally ascending (m a+ i = m a + l), flat (m a+ i = 
m a ), or descending (m a+ i = m a — 1): 




The resulting quiver encodes the skew-symmetric matrix B m . 

Finally, let us mention the following Lemma (see Ref.[3] for details), used crucially in 
the following. 

Lemma 2.1. Any Motzkin path m with rrti > for all i £ [l,r] may be attained from mo 
via iteration of forward mutations of the form fif acting at each intermediate step on a 
Motzkin path m in either of the two following local configurations around i: 

• Case (i): rrii-i = mi = m i+ i — 1 

• Case (ii): mj_i = rrij = m i+ i 

2.4.3. Quantum cluster algebra for the quantum A r Q-system. The quantum cluster al- 
gebra corresponding to our quantum Q-system (12.1 1) has the fundamental seed x = 
((Ri,o)i=i, (Ri,i)i=i), and the same exchange matrix B as in the commuting case. The 
commutation relations (12. 2p correspond to taking the initial compatible pair (Ao,I?o) such 
that A = (r + l)!?^ 1 . The compatibility implies that A m = (r + 1)!?^ 1 for all Motzkin paths 
m, and we have explicitly (AJ^i^j = {jp — m)Ajj in terms of the matrix A of eq. fl2.3j) . 
for all pairs (i, m) and (j,p) of cluster indices in x m , which leads to the commutations (12. 2p . 

3. Quantum system solution for R 1<n via continued fractions 

3.1. Generating functions. We set R n = i?i iTl . To each Motzkin path m = (m 1; m 2 , m r ), 
we associate the generating function 

oo 

(3-1) F m (t) = y^ j t n R n+mi R^ n \ 

n=0 

We also use the "rerooted" generating function 

oo 

(3-2) G m (t) = 2^ t n Rn+mx+lRrrn+l 

n=0 

3.2. Continued fraction expressions. For some variables x±, X2, X2 r +i elements of 
a non-commuting algebra, let Ji(xi, X2, X2 r +i) be the "non- commutative Jacobi-type 
(finite) continued fraction" defined inductively by 

Js{%2s-l, %2s, ■1 %2r+l) = (l — x 2s-l ~ Js+1 (^2s+l , ^2s+2 j • • • j ^2r+l ) x 2s) (1 < S < r + 1) 

Jr+2 = 
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112 13 r 1 r+1 

FIGURE 1. The weighted graph T r , with r + 1 vertices labeled 1,2, ...,r + 1. 
For each oriented edge e, we have indicated the corresponding step weight w(e), 
namely: w{i —>•«) = X2i-\ (i = 1, 2, r + 1), 10 (i — > i + 1) = X2i, w(i + 1 —>•«) = 1 
(i = l,2,...,r). 



For some central scalar parameter i, the function J\(txi, ...,tx2 r +i), expanded as a formal 
power series of t, may be interpreted combinatorially as the generating series for "quantum" 
paths on the weighted graph T r depicted in FigfTJ from and to the origin vertex 1, each 
path being weighted by the product of its step weights taken in the order they are traversed 
and with an extra weight t per step along each edge pointing toward the origin. 

Definition 3.1. To each Motzkin path m we attach a sequence of "weights" y(m) = 
(yi (m), ?/2r+i(rci)) via the following induction under forward mutation fj,f(va) = m' , 
depending on whether m is in cases (i) or (ii) of Lemma \2. 1\ above. For short we write 
yi(m) = yi and yi{va!) = y[. First, we have y'j = yj for j ^ 2i — 1, 2i, 2i + 1 (case (%)) and 
j 7^ 2i — 1, 2i, 2i + 1, 2i + 2 (case (ii)), while: 

{y-H-i = UK-* + 
y r 2i = y2i+mi(y2i-i + y2i)~ l 
y'21+1 = y2i+iV2t-i{y2t-i + y2i)~ 1 

(3.3) Case(ii) : y' 2l+2 = j/m+iJ/m-iCj/m-i + 2/2i) _1 

This determines the y's entirely in terms of the initial data y(m ). We have the following 
Theorems. 

Theorem 3.2. (^7\) For the fundamental initial data with m = m = (0,0,...,0) ; the 
solution to the quantum A r Q-system satisfies the following identity: 

(3.4) F mo (t) = 1 +tG mo (t)y 1 (va ), G mo (t) = Ji(%(m ), ty 2 r+i{m )) 
where 

(3.5) 2/2i-i(m ) = Ri^R^^RilRi-ifl y 2 i(nio) = Ri+i^R'J R^q Ri-i^ 
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Theorem 3.3. (^\) For y(m ) = (yi(m ), ...,y 2r+ i(m )) given by (13.51) . the solution of 
the recursion (I3.3P reads: 

(3.6) 2/2i-i(m) = g*" 1 i?i imi+ ii? l "^i?i_i imi _ 1 i? i l 1 1/m ._ i+1 

(m) = /^ +1 ' m *+i +1 ^+ 1 i"»i+i^+ 1 > m i^+ 1 .^i+i ^ m * = m * +1 ~^ ^ 
I 1 otherwise 



(3.7) 



x Ri+l,mi + lRi jmi +lRi trni Ri-l,mi 

^i-l,m i +lRi-l t m i Ri-l,m i ^ 1 Ri-l t m i _ 1 +l if m i = ~~ 1 

1 otherwise 



Theorem 3.4. (]7]j For £/ie |/j(m) as in Theorem VJ.^ we have: 

F m (t) = 1 + tG m (t) yi (m), G m (t) = J x (y i(m), y 2r+1 (m)) 
where % = ?/j(m) are related to yi = yi(m) via: 

{Vn-i = ty 2i -i y 2i = ty 2i if m i+1 = mi 

2/2t-i = t(y 2i -i + y 2i ) y 2i = t 2 y 2i+1 y 2i if m i+1 = m* + 1 (ie [1, 2r + 1]) 
jtet-i = ty-u-i - y 2l l +1 y 2 i m = y 2 l+\yu if m+i = m - 1 

Example 3.5. For m = m = (0, 0, 0), we have % = tyi for i = 1, 2, 2r + 1. The 

generating function reads 

G mo (t) = (1 - % - t(l - % - t(...(l - fy 2r _! - t(l - ty 2r+l y 1 y 2r )- 1 ...)- 1 y A )- 1 y 2 y 1 

Example 3.6. For m = mi = (0, 1, ...,r — 1), we have y 2 %-\ = t(y 2 i-i + y 2 i) and y 2 i = 
t 2 y 2 i+iy 2 i for all i, where yi = i/j(mi). The result is simplest for F mi (t) = 1 + tG mi (t)yi. 
The generating function may be rearranged using the following identity at each step: 

(3.9) a + 6 + (1-c - u)~ 1 cb = a + (1 - c - u)~ l {l - u)b = a + (1 - (1 - u) _1 c) -1 6 

The result reads 

F mi (t) = (1-*(1- t(...(l - £(1 - tSHil'V)' 1 -)' 1 *)"^)" 1 

27ms zs easi/?/ expressed in terms of the "non-commutative Stieltjes-type (finite) continued 
fraction" defined inductively as 

S k (x k , x k +i, x 2r+ i) = (1 - S k+1 (x k+h x 2r+ i)x k y l (1 < k < 2r + 1), SW+2 = 1 
w'a: F mi (t) = Si(ty 1 ,ty 2 , ...,ty 2r+1 ). 

In view of the Example 13 .61 we may write another (mixed Stieltjes-Jacobi-type) continued 
fraction expression for G m (t) for arbitrary m. This will be crucially used in the following. 

Any Motzkin path m = (m 1; m r ) may be decomposed into strictly ascending segments 
of the form (m,m + 1, ...,m + k — 1) separated by weakly descending steps of the form 
(m, m) or (m,m — 1). Accordingly, we may transform the Jacobi-type fraction for G m (t) 
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in Theorem 13.41 by "undoing" the pieces of the continued fraction that correspond to the 
strictly ascending segments of m. 

To this end, assume that we have a strictly ascending segment (m io ,m io+1 , ...,m io+k _i) 
with m io+ j = m io + j for j = 0, 1, k — 1, which is followed by a weakly decreasing step 
{m io+k ^m io+k ) } with m io+k = m io+k - X or m io+k = m io+fc _i - 1. Then by definition, we 
have the relations 

Ji +j — (l ~~ t(y2i +2j-l + U2i +2j) — t 2 Ji 0+ j + iy 2 i +2j+iy2io+2j) (j £ [0, k — 2]) 

(3.10) J io+ k-l = (l - ty 2 i +2k-3 - (Ji +k - $m t0+k ,m l0+k _ 1 -l)y2i +2k-2) 

Like in Example 13.51 let us rearrange this piece of fraction iteratively from the bottom up, 
by using the relation (13. 9p . This allows to rewrite: 



Jio — (1 ~ ty2i -i — tS 2 i y2io) 1 
S 2io +j = (1 - tS 2io+j+ iy2i 0+ i +j ) _1 (j e [0, 2fc - 4]) 

£2i +2fc-3 = (l — S 2 i +2k-2y2i +2k-2) 
(3.11) 5'2i 0+ 2fc-2 = ^i +fe - ( ^m io+fc ,m io+fc _ 1 -l 

Indeed, we start by using (13.91) at the step j = k — 2 of (I3.10p . with a = ty 2 i + 2 k-5, 
b = ty 2 i 0+2 k-4, c = y 2 i +2k-3 and u = S 2io+2k - 2 y 2iQ+2k - 2 , to rewrite: 



Ji +k-2 — (l — ty 2 i Q + 2k -5 — t(l — t(l — S 2 i +2k-2y2i +2k-2) 1 y2i +2k-s) 2/2 



1 

i +2fc-4 



Iterating this leads straightforwardly to (13. lip . Repeating this transformation for every 
strictly ascending segment of m, we arrive at: 

Theorem 3.7. For any Motzkin path m ; we have the following mixed Stieltjes-Jacobi type 
continued fraction expression for G m . Leti±, ..,i s andii, ...,£ s be the positions and lengths of 
the strictly ascending segments of m, of the form (m ia ,m ia + 1, ...,m ia+ £ a _ 1 ) , a — 1,2, ...,s. 
Then we have: 

G m (t) = J\ 

J%= (1 - fai-i ~ Ji+iV2i) 1 (« £ U a [i tt , i a +4 - 1]) 

■4 = (l - ty 2ia -i - S 2i ty2i a y l (a G [1, s\) 

Si= (1 - S i+l ty i+i y l (i G [2i a ,2i a + 24-4];ae [l,s]) 

5*2^+2^-3 = (l — £2^+2^-22/2^+2^-2) (ae[l,s]) 

5 , 2i a +2£ a -2 = ^i a +^ a - ^m ia+ « a ,mi a+£a _ 1 -l (a G [1, s]) 



10 PHILIPPE DI FRANCESCO 

3.3. Quantum commutation relations for the weights. Using the commutations 
( 12. 2p . we obtain: 

Theorem 3.8. Introducing p = q r+1 , the weights y^ = yj(m ) of Theorem \ 3.S\ obey the 
following p-commutation relations: 

ViVi+i = PVi+iVi (i e [l,2r]) 

V2iV2i+2 = py2i+2V2i (ie[i,r-i]) 

= otherwise 

Using the recursion relations (13. 3p . we deduce the commutation relations of yi(m), for 
arbitrary Motzkin paths m: 

Theorem 3.9. For a given Motzkin path m, the weights yi = y%{va) have the following 
commutation relations 

(3.12) yiy i+1 = py i+ iyi (ie[l,2r]) 

(3.13) y 2i y 2i+2 = y 2i+2 y 2i (ie[l,r-l]) 

VjVk = Vkyj otherwise 

Proof. By induction under mutation. The Theorem holds for m = mo (Theorem 13. 8j) . 
Assume it holds for some m, then consider m' = pi (m) in either cases (i) or (ii) of Lemma 
12. 1[ and denote by y'j = ^ (m'). We deduce the following commutations from the recursion 
hypothesis (eq. (13.121) ): 

(3.14) y 2 j-iy 2 k-i = V2k-iV2j-i (j, k G [1, r + 1]) 

(3.15) (y 2j -i + y2j)y2j+iV2j = py 2j+ iy 2j (y 2j -i + y 2j ) (je[l,r\) 

(3.16) {y2j-i + y2j)y2j+iV2j-i = y 2j+ iy 2j -i{y 2j -i +y 2j ) (je[l,r)) 

In both cases (i) and (ii), this implies that y' 2 j-\y 2 k-i — v'lk-xV^j-x ^ or a ^ 3ik and that, 
as rm-i = rm i n bot h cases, j4_ 2 ^-i = PVn-iV^i^ w hile y '2i -\y'2i = PV-hVm-i as a 
consequence of (13.151) and y 2i y 2i+ i = P2/2i+iZ/2i by use °f (I3.15p - 03.16p . Finally, we get 
y'2i-2V'2i = y'2iV'2i-2 as y'21-2 = V2i-2 p-commutes with both y 2i _ x + y 2i and y 2i . Using the 
expressions (I3.3P for the cases (i) and (ii), we finally find 

• Case (i): y' 2i y' 2i+2 = py 2i+2 y 2i , as y' 2i+2 = y 2i+2 and y 2 iy2i+2 = 2/2i+22/2i (horn m i+1 = 
mi + 1). 

• Case (ii): y' 2i y' 2i+2 = P 2 y 2i+2 y 2i , as (y 2 i-i + Uk) commutes with y 2i + 2 y 2i -x (due to 
2/2i2/2i+2 =py2i+2y2i, from m i+ i = rrii). 

Noting finally that m' i _ l — m • + 1 = in both cases (i)-(ii), and that m! i — m' i+l + 1 = 1 in 
case (i) and = 2 in case (ii), the Theorem follows. □ 

To each Motzkin path m we may associate a quiver Q m with 2r + 1 vertices labelled 
i = 1,2, ...,2r + 1, that summarizes the commutations of the ?/j(m)'s as follows: we draw 
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FIGURE 2. The nine quivers Q m corresponding to the nine Motzkin paths m = 
(mi,m2,m3) G M3, and Q mo =(o,o,o) 011 top of the diagram. In each quiver, the 
vertices are all labeled from bottom to top 1, 2, 7. The arrows between quivers, 
labeled i, correspond to mutations fif in either case (i) or (ii). 



m arrows from vertex i to vertex j whenever yiUj = p m yji/i. For illustration of Theorem 
13. 9\ we have depicted in FigfS] the example r = 3 of the quivers Q m in the fundamental 
domain M3 of Motzkin paths under global integer translations, and indicated by arrows 
and superscripts i the mutations \i\ acting on them. 



3.4. Quantum multinomial expressions. 



3.4.1. p- combinatorics. For a\, ...,ctfc £ the quantum multinomial coefficient is denned 
as: 

j-jaiH ha fc j 



ai + 
ax, 



■ + a k 



Li=l 



UUUZxii-P 1 
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We assume by convention that when k = 2 







1. We have the following p- 



p 



multinomial identity for k variables xi, such that XiXj = pxjXi for all 1 < i < j < k: 



(xi H h x k ) 



E 



m-i ,...,,m^,>0 
m| H mp.—n 



n 

m 1 ,...,m k 



x k x k-l x l 



This reduces to the standard p-binomial identity for k = 2. 
We also define the following formal generating series: 



i=0 



m=0 



+ m — 1 
m 



and, for variables Xi, such that x^Xj = pxjx^ for all 1 < i < j < k, we have: 



(x 1 + ... + x k )= J2 

mi,...mfe>0 



1 + mi H \-m k 

'-l,m 1 ,...,m k 



• x 1 



X k X k-1 



3.4.2. Explicit expressions for G m (t) . We have the following explicit expression for G mo (t), 
after ordering of the yi = yj(m )'s. 

Theorem 3.10. For the flat Motzkin path mo = (0, 0, 0), we have: 



G mo (t) 



e n 

«l,* 2 .-".*2r+l>° i=0 
*0=V2r+2=0 



— 1, @2i+l, ^2i+2 



(ty 2r+i y^(ty 2r )^...(t yi ) 



1 + 



e n 

<l>l,i S ,...,i2r+l>0 1=0 



%i + ^2i+l + ^2i+2 — 1 
C-2i — 1, ^2i+l) ^2i+2 



(^ 2r+1 )^ i (^) fe, ----(%r i 



*0 = M 2 r+2=0 



Proof. By induction. We start with the expression (13.41) for G mo (t) = Ji, and write the 
inductive definition J x = (1 —ty 1 — t^ite) -1 - Next we note that y^JzV?) = vi^yzlVi-, as J2 
only involves i > 3, which all commute with j/!, by Theorem 13 .81 We deduce the formal 
expansion 

£2 



^2>Q 



p 



Assume we have 

(3.i7) Ji= e n 



TO— 1 r 



<l.->*2m>° i=0 
«0=1 



?2i + ^+1 + ki+2 - 1 
@2i — 1, ^2i+l, ^2i+2 
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for some m > 1, then by Theorem ESI we know that ?/2m2/2m+i = M/2m+i?/2m, Z/2m?/2m+2 = 
Py2m+2y2m while y 2 m commutes with y i5 i > 2m + 3, and in particular with J m +2- This 
implies: 

y2m-Jm+l = 2/2m(l "~ ty2m+l — tJ m +2y2m+2) 1 = (1 ~ Pty2m+1 ~ ptJ m +2y2m+2) 2/2m 



r. 



<j>e(ty2m+l + tJ m+ 2y2m+2){ty2m) 

^ + 4m+l + 4m+2 — 1 
- — 1, ^2m+l) 4m+2 



E 

^2m+l,^2m + 2>0 



(t Jm+2y2m+2Y 2m+2 (ty2m+lY 2m+1 (ty2rn 



where in the last line we have used the fact that ?/2m+i2/2m+2 = py 2 m+ 2 y 2 m+i> while ?/2m+i 
commutes with J m+ 2- Substituting this into (I3.17p . with £ = £ 2 m, we obtain the same 
summation with m — >■ m + 1. We conclude that (13 . 1 7[) holds for all m, and in particular for 
m — r + 1, in which case J r+ 2 = imposes that the last summation reduce to 4r+2 = 0, 
and the first part of the Theorem follows. The second follows trivially from the relation 
between F mo and G mo . □ 

Theorem 3.11. For the ascending Motzkin path mi = (0, 1, 2, r — 1), we have: 

2r <~i i) 1 r 1 

(ty 2r+1 ) e ^(ty 2r )^...(t yi yi 



G mi (t) 

Fm,(t) 



e n 

«l,*2,.»rfer+l>0 i=l 
2r 

e n 

fi,fc,..,&+i>0i=l 



?i + 4+1 - 1 + <fi,l 
4+1 



+ 4+1-1 
4+1 



(^2r+l)^ +1 (^)^---(%) 4 



Proof. We start from the expression F mi (t) = Si of Example 13.6} and write 

(^3 2 i/2)(%; 



Si = (1 - ts 2 yi)- 1 = MSm){tyi) h = J2 

4>0 <1^2>0 



? i + 4-i 

4 



Assume we have 

*= e n 



m-2 



tl,...,£ m _l 



>0 i=l 



+ 4+1 
4+i 



(tS m y 



m—l i 



-'(ty 



m-2 



■ (tyxY 1 



for some m, then writing {tS m y m -iY = 0^(tS' m+ i?/ m )?/^„ 1 immediately implies the same 
relation for m — > m + 1. We conclude that it holds for all m, in particular for m = r + 2, 
where it boils down to the second part of the Theorem. The first part follows from the 
relation between F mi and G mi . □ 



3.4.3. The main Theorem. More generally, we have 



14 PHILIPPE DI FRANCESCO 

Theorem 3.12. For a generic Motzkin path m, we have the following: 
F m (t) = 



G m (t) 



J2 A m (£ u £ 2 , ...J 2r+1 )(ty 2r+1 ) e ^ . . . ( tyi ) h 
i 1 ,e 2 ,...,e 2 r+i&+ 

^m(il + M 2 , ...,£2r+l){ty2r+lY 2r+1 " " " (tyl) 



where A m is defined as the product 
A m (£i, ...,£ 2 k+i) 1 



i: 



r-2 

n 

i=l 



%% + ^2i+l + ^2i+2 — 1 
^2i — 1) ^2i+l) ^2i+2 



if m i+1 = mi 



2i+1 



? 2i + ^2i+l + hi+2 ~ 1 
^-1 



hi + £2i+l - 1 
^2i+l 



^2i + ^2i+l + ^2i+2 
^2i+2 



?2i+l + ^2i+2 — 1 
^2i+2 



if irij 



i+1 



rrii - 1 



2i+1 



if mj+i = mj + 1 



2i+1 



-2r + ^2r+l — 1 
^2r+l 



• 2r+1 
1 2r 



2i+2 ^ 




where we have represented the local structure of the corresponding quiver Q m that encodes 
the p- commutations of the y 's. 

Proof. As before, we proceed by descending induction. Assume that for some k > 1, such 
that rrik = m k -i (Case (1)) or m k = rrik-i — 1 (Case (2)), we have an expression of the 
form: 

F m (t) = Yl A 

mi,...,m k ■■■} <-2k) 

h,...,l 2k &% + 

x ((J fc+ i - e k )y 2k ) t2k (ty 2k ^y^ ■ ■ ■ 

£k — r , 



(3.18) 
(3.19) 

We now have three possibilities for the Motzkin path: 
Case (a): m^+i = m k . 



J m k _ 1 ,m k +l 
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Case (b): m k+ i = m k — 1. 

Case (c): m k+ j = m k + j for j — 0, 1, 2, £ — 1 and < m^_i. 

Writing J fe+1 = (1 - y 2k+1 - 4+2^+2) _1 , we note that we always have y 2 kV2k+i = 
pfak+ihk, V2kV2k+2 = PV2k+2V2k, while y 2k commutes with J k+2 , as the latter only depends 
on yi with i > 2k + 3. This allows to write: 

Case(l): (J k +iy2k) m = <j> m (fak+i + ^+2^+2) (^/2fc) m 

Case(2): ((J fc +i - l)?/2fc) m = p m( ™ 1 m (y 2fc +i + 4+2^+2)^+1 + 4+2^+2) m (2/2fc) m 

(3.20) = <fr m (y2k+l + Jk+2y2k+2){y2k+l + J ^+2^+2)"' '(*2/2fc+l)~ W '(^fc)™ 

where we have used y 2 fe = ^2/2fc m the case (1), and ^ 2 fc = 2/2fe+i2/2fc in the case (2), hence 

fer = ^ I ^(^2 fc +i)~ m (^) m - 

We first treat the Cases (a) and (b). In both cases, we have y 2 k+W2k+2 = Pll2k+2y2k+i, 
while ?/2fc+i commutes with Jfc +2 . This suggests to rewrite 

2/ 2 fc+l + Jk+2y2k+2 = ty 2 k+l + (Jk+2 — eA+l)?/2A+2, Cfc+1 = ^m fc+1 ,m fc -l 

in which the two summands p-commute. In the Cases (la) and (lb) we get: 



4>m(ty 2k +l + (Jk+2 ~ tk)y2k+2) = ^ 

j,n>0 

while in the Cases (2a) and (2b): 



m — 1 + j + n 
m — n 



((Jk+2 - £k)y2k+2) n (ty 2 k+iy 



4>m(ty2k+l + (Jk+2 - ^k)hk+2)(ty2k+l + (Jk+2 ~ ^k)y2k+2) 

m — 1 + j 

j 



E 

j>0 



(ty2k+l + (Jk+2 — e k)y2k+2 



,m+j 



E 

j,n>Q 



j + n — 1 




j + n 


m — 1 







((4+2 - £k)y2k+2) n (ty 2 k+iy 



We now deal with the Case (c). By Theorem 13. 7\ the strictly ascending segment of length 
starting at m k corresponds to the mixed Stieltjes-Jacobi expression: 



Jk+l — (1 — ty 2 k+l — S 2 k+ 2 ty2k+2) 



1-1 



Si — (1 — Si + ityi + i r 
S 2 k+2i-i = (1 — (Jk+e+i — £k+e)y2k+2e, 



(i e [2k + 2,2k + 2£ - 2]) 

-l 
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Eq.f l3.20p may be rephrased in the present case by substituting y2k+i + Jk+2V2k with ty2k+r 
S 2 k+2ty2k+2- As before the two cases (1) and (2) lead respectively to: 



(j>m{ty2k+l + {S2k+2)ty2k+2) = 

j,n>0 L 

4>m{ty2k+l + {S 2 k+2)ty2k+2){ty2k+l + {S 2 k+2)ty2k+2) 



m — 1 + j + n 
m — 1, j, n 



{S2k+2ty2k+2) n {ty2k+i)' J 



V 
rn 



E 

j,n>0 
j-\-n>m 



j + n — 1 




j + n 


m — 1 


p 


j 



(S2k+2ty2k+2) n (ty2k+i) J 



Due to the commutation relations between the y's involved, we may apply directly the 
recursion in the proof of Theorem 13.111 to get: 



(5 , 2fc+2^2/2fc+2j 



E 



2fc+2£-l r 



n 



«2fc+3' £ 2fc+4>---^2fc+2£>° i=2fc+2 
e 2k+2= n 

x ((Jk+e+i — £k+e)y2k+2£ / 



+4+1-1 
4+1 



t2k+2t (ty2k+2£-xY 2k+2e - 1 ---(ty: 



2k+3, 



t2k+3 {tV2k + 2) 



Assembling all the cases, we find the following recursion relation for A in the cases (ia) 
and (ib), with i = 1,2: 



(3.21) 



Am\,...,m k +i (4; ■••■> 4ft+2) y 4mi,...,m fc 4fe+l, 4fc+2j 

2k + 4fc+l + 4fc+2 _ 1 
4fc — 1) 4fc+l> 4*+ 2 

2k + 4fc+l + 4fc+2 _ 1 
£2k~l 



u. 



j p 

() -2k + 4fc+l + 4fc+2 
4fc+2 



and in the cases (ic), i = 1,2: 
(3.22) .4 



mi,..,m t+w (4, ■■•,4fc+2^) 

4&+1) 4*+2)^(4a;+2, 4fc+2l) 

2fc+2^-l 



+ 



-i+1 



1 



^+1 



(3.23) V r (£ 2fc+2 ,...,£ 2fc+2 ,)= 

i=2£+5 

This determines A entirely, with the initial conditions: 

Case (a) or (b) : A mi (4, 4 



i 



Case(c) : A mi 



m2,...,m< 



<«>= n 



i=l 



+ 4+1 ~ 1 
4+1 



P 
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and the Theorem follows, with A m (£i, izr+i) — ^mi,...,mr,m r (4, • ^2r+i, 0), as the last 
step has ?/2r+2 = 0. □ 

Extracting the coefficient of t n in the series of Theorem 13.121 and using the expressions 
for the weights given by Theorem 13.31 we immediately deduce: 

Corollary 3.13. The solution n of the quantum A r Q-system is expressed as a Laurent 
polynomial of any admissible initial data, with coefficients in Z + [g,g -1 ]. 
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